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Abstract 

Satellite-based  navigation  systems,  such  as  GPS,  Glonass,  and  the  planned  Galileo 
European  system  rely  on  stable  atomic  clocks  for  their  operation.  Space-borne  and  ground- 
based  clocks  provide  the  time  frame  allowing  the  users  to  perform  passive  ranging 
measurements  and  quickly  locate  their  position  [13]. 

From  the  operational  point  of  view,  timing  stability  enters  the  overall  accuracy  budget 
together  with  orbit  estimation  and  operational  constraints.  User  ranging  accuracy  is  a  complex 
function  of  many  factors,  and  optimum  selection  of  the  clocks  requirements  is  of  paramount 
importance  to  optimize  the  overall  system.  Ideally,  user  range  accuracy  is  affected  by  both 
timing  and  orbit  estimation  accuracy.  In  turn,  both  timing  and  orbit  estimation  are  affected  by 
the  oscillator  intrinsic  accuracy  and  stability.  Stability  plays  a  role  first  in  determining  the 
optimum  time  interval  over  which  the  estimation  is  carried  on  and,  secondly,  in  determining  the 
maximum  length  of  the  time  interval  over  which  the  computed  estimates  can  be  propagated  and 
still  be  within  the  navigation  system  specifications. 

Operational  constraints,  such  as  the  maximum  uploading  rate  for  civilian  systems,  and  the 
maximum  unattended  system  operation  vs.  degraded  performance  for  military  systems,  interact 
with  purely  technical  considerations  to  define  the  “best”  clock  for  a  given  application.  To  all 
this,  we  should  add  the  requirements  placed  on  the  ground  clocks  which  play  a  significant  role 
in  determining  the  overall  system  performance  especially  in  systems  (such  as  GPS  or  the 
planned  Galileo)  relying  on  the  one-way  passive  ranging  technique  for  orbit  and  clocks 
maintenance. 

In  the  paper,  we  will  discuss  the  requirements  forced  on  the  clocks  by  the  specific 
application,  as  well  as  the  consequent  tradeoff  dictated  by  the  stability  of  the  orbit  and 
operational  constraints.  The  clock  constraints  dictated  by  the  simultaneous  requirement  to 
produce  and  disseminate  timing  information  related  to  UTC  (or  TAI)  will  not  be  discussed,  as 
we  will  limit  our  discussion  to  the  aspects  related  to  the  navigation  services  only.  In  this  frame, 
we  would  like  also  to  present  the  rationale  for  the  engineering  development  of  a  novel  frequency 
standard  (a  Rb  passive  maser)  being  developed  by  the  IEN  and  Alenia  Spazio  under  an  ASI 
(Agenzia  Spaziale  Italiana)  contract  specifically  for  possible  applications  in  future  navigation 
systems. 

The  opinions  expressed  in  the  papers  are  of  the  authors  and  do  not  imply  necessarily  the 
endorsement  of  the  respective  organizations. 
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SYSTEM  MAINTENANCE 


Satellite-based  navigation  systems  use  one-way  ranging  measurements  for  system  orbit  estimation  and 
time-keeping.  The  operational  advantage  of  one-way  ranging  versus  two-way  ranging  is  obvious  when  one 
considers  the  complexity  of  the  ground  monitoring  stations  in  the  two  approaches.  One-way  requires  a 
simple  omni-directional  antenna  to  track  all  the  satellites  in  view,  is  completely  passive  (non-transmitting) 
and  the  station  can  be  deployed  or  re-deployed  with  minimum  effort,  requiring  only  a  surveyed  location, 
making  it  ideal  for  a  military  system.  On  the  contrary,  a  two-way  ranging  station  requires  a  complex 
transmitting  equipment,  a  large  directional  antenna  and,  as  a  consequence,  will  not  be  able  to 
simultaneously  track  multiple  satellites,  and  it  is  expensive  to  deploy. 

One-way  ranging  measurements  are  termed  “pseudo-ranges,”  since  they  contain  the  system  clocks 
contributions,  namely  the  space  vehicle  (SV)  clock  and  the  monitor  station  (MS)  clock  in  addition  to  the 
propagation  delay  caused  by  finite  propagation  velocity  vp  over  the  range: 

(1)  P  =  Vp  '  (Af prop(SV-MS)  +  A* ms-SV  )=Vp  ‘  I.A t prop(SV-MS)  +  fasv  +  ^  MS  )J 

where  in  the  last  term  we  have  expressed  the  contribution  A tMS_sv  to  the  pseudo-range  measurement  as 

the  sum  of  the  offsets  of  the  individual  clocks  with  respect  to  the  system  time  which,  for  a  composite  clock 
solution1,  must  satisfy,  in  principle,  the  relationship: 

(2)  ^  '  Af sv,,  +  ^  Wj  •  A tMS  j  =  0 

i  j 

and  the  sums  are  carried  on  over  all  SVs  and  MSs  clocks  in  the  system,  each  properly  weighted  with 
weights  w,  and  wr 

The  fundamental  assumption  of  the  one-way  ranging  technique  is  the  capability  to  separate  the  three 
contributions  to  the  pseudorange  measurements,  given  a  sufficient  number  of  measurements.  This  will 
yield  range  observables,  used  to  update  the  estimate  of  the  orbit,  and  time  offsets  observables,  to  estimate 
the  clock  offsets  and  derive  from  the  latter  the  clock  parameters  (phase  and  frequency  offsets  and  drift) 
subject  to  condition  (2). 

However,  a  sufficient  number  of  independent  measurements  are  required  to  solve  the  system  unknowns, 
but  additional  constraints  apply  for  the  term  separation  to  be  effective,  namely  that  a  sufficient  number  of 
simultaneous  independent  measurements  be  available.  Simultaneity  of  observations  from  the  same 
monitor  stations  of  a  large  number  of  satellites,  as  well  as  simultaneous  visibility  of  multiple  stations  by  an 
individual  satellite,  is  the  key  to  an  effective  separation  in  the  solution  of  the  clock  contributions  from  the 
pseudo-ranges. 

Single  and  double-differencing  techniques,  widely  used  for  data  reduction  in  the  geodetic  community,  will 
be  of  help  in  understanding  the  underlying  physical  rationale.  When  two  satellites  are  simultaneously 
observed  by  a  single  monitor  station  (Figure  1),  the  first  difference  of  the  two  pseudoranges  drops  the 
common  MS  clock  term  from  the  observable. 


1  If  a  Master  Clock  solution  is  chosen,  all  time  offsets  will  be  referenced  to  this  single  clock  (or  ensemble). 
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Figure  i  -  Simultaneous  observations  from  a  single  monitor  station 


Notice  that  we  assume  that  the  ionospheric  propagation  effects  are  completely  removed  by  the  use  of  the 
dual-frequency  technique  and  tropospheric  effects  are  common  to  the  two  measurements",  so  they  cancel 
out  too.  This  yields  a  nice  observable  for  the  SV  clocks.  When  a  single  satellite  is  simultaneous  in  view  of 
two  ground  stations,  the  situation  depicted  in  Figure  2  applies. 


Figure  2  -  Single  satellite  in  simultaneous  view  of  two  monitor  stations 


Again,  the  first  difference  drops  the  SV  clock  and  yields  an  observable  which  contains  only  the  MS  clock's 
contribution.  However,  while  the  same  considerations  for  the  ionospheric  propagation  still  applies  as 
before,  now  the  tropospheric  delay  is  not  "common  mode"  and  will  affect  the  final  estimation  of  the  range 
and  MS  clock  estimation. 

Monitor  stations  clocks  estimation  (prediction)  is  affected  by  ’local"  propagation  delays  which  are  not 
correlated.  Hence,  errors  propagate  to  orbit  determination  and  indirectly  affect  the  final  user 
positioning/timing  accuracy. 

Since  a  second  difference  of  a  number  of  simultaneous  observations  having  in  common  the  SVs  and  the 
MSs  will  yield  the  orbit  estimation  free  of  clock  terms,  it  is  intuitive  that  the  tropospheric  effects  are  the 


'  This  may  not  be  always  the  case,  since  local  perturbations  to  the  atmospheric  component  of  water  vapor  and 
temperature  along  the  line  of  sight  lend  to  be  unequal  when  looking  in  different  directions.  In  this  respect,  there 
are  contrasting  opinions  on  the  requirements  for  additional  measurements  along  the  line  of  sight  or  concerning  the 
magnitude  of  the  effect  with  respect  to  modeling;  but  certainly  the  change  in  the  index  of  refraction  ol  the 
troposphere  can  be  major  when  considering  geographically  separated  sites. 
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major  source  of  error  left  in  the  MS  clocks  state  estimate  and,  as  a  consequence,  of  the  MS  clock  states 
error  projection  on  the  orbit  estimate. 

Therefore,  use  and  improvement  of  meteorological  data  and  tropospheric  propagation  models  is  of 
importance  in  the  overall  system  error  budget.  In  parallel,  an  independent  monitoring  capability  of  the  MS 
clock  behavior  (by  two-way  time  transfer,  for  instance,  as  if  has  been  proposed  in  several  occasions  in  the 
past  [14])  may  help  in  highlighting  possible  mismodeling  effects  in  the  troposphere,  as  well  as  improving 
the  capability  to  verify  '  the  MS  clock'  s  state  estimate  and  its  final  contribution  toward  the  orbit  estimation 
errors* * * 4. 


Looking  at  the  overall  system 
data  flow,  some  considerations 
are  readily  apparent  (Figure  3): 
measurements  errors  (pseudo¬ 
ranges)  contains  SV  and  MS 
clocks  terms,  as  well  as  range 
(orbit-dependent)  terms.  The 
measurement  errors  are  mainly 
radial  and  the  OD&TS^  process 
projects  these  errors  into  the 
various  components  of  the 
system  state  vector  to  yield 
estimates  of  SV  clocks,  MS 
clocks,  and  orbit  elements,  all 
of  which  are  affected  by  the 
projections  on  the  various 
states  of  the  essentially  radial 
measurement  errors. 
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Figure  3  -  Overall  navigation  system  data  flow 


User  navigation  solution 


USER  POSITIONING 
SOLUTION 


The  navigation  solution  in  the  user  set  recombines  SV  clocks  and  range  estimates  based  on  system  state 
prediction,  where  all  errors  recombine  again  in  the  radial  direction.  Since  the  user  clock  behavior  can  be 
neglected  when  solving  for  the  PVT  (position,  velocity,  and  lime)  user  vector,  the  strong  correlation  of  SV 
and  orbit  estimation  errors  in  the  radio  I  direction  is  beneficial  for  the  user.  However.  MS  stale  estimate 


As  an  experimental  device,  a  very  useful  technique  involves  the  use  of  LASER  ranging  to  a  retrore  flee  tor  array 
placed  on  the  satellite.  Besides  the  obvious  result  of  verification  of  the  ranging  measurements,  it  is  to  be  recalled 
that  the  troposphere  is  a  dispersive  medium  at  optical  frequencies,  so  that  a  two-color  ranging  system  may  yield 
information  about  the  troposphere  refraction  properties.  The  resulting  data  can  be  utilized  also  to  correct  the  radio 
frequency  measurements  along  the  line-of-sight  to  the  SV.  The  final  addition  of  a  two-color  LASSO-type  detector 

onboard,  which  is  adding  very  little  in  terms  of  additional  weight  and  power  consumption,  will  allow  the  full 

separation  of  all  the  relevant  terms  involved  in  Eq.  (1)  with  higher  accuracy  that  available  from  the  current 
operational  measurement  techniques. 

4  As  a  possible  alternative,  two-way  data  may  contribute  to  the  frequency  offset  and  drift  estimation  ot  the  MS 
clocks,  leaving  the  time  offset  determination  to  the  GDT&TS  function  alone  (to  satisfy  Eqs.  (1)  and  (2) 
simultaneously). 

OD<&TS  -  Orbit  Determination  and  Time  Synchronization  function. 
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errors  do  not  enter  the  user  solution;  therefore,  their  effect  on  the  estimation  of  the  SVs  clocks  and  orbit 
will  affect  the  final  user  prediction  embedded  in  the  PVT  solution. 

The  previous  considerations  lead  to  the  use  of  stable  oscillators  in  the  ground  monitor  stations,  where 
stability  is  to  be  emphasized  in  the  double  areas  of: 

•  short-  to  medium-term  stability,  not  to  degrade  the  measurements  over  the  time  period  used  for  orbit 
updating,  as  well  as 

•  long-term  stability  to  allow  a  clean  separation  of  the  clocks  states,  with  the  aim  of  improving  the 
navigation  system  time  stability  and,  overall,  its  orbital  estimation  accuracy. 


ORBIT  MAINTENANCE 

To  maintain  the  navigation  accuracy  of  the  system,  orbit  and  time  parameters  need  to  be  periodically 
updated,  since  both  tend  to  drift  with  time  and  to  increase  the  related  errors.  The  periodical  maintenance  is 
based  on  frequent  updates  of  the  navigation  parameters  and  less  frequent  tuning  of  system  parameters  or 
hardware  calibration  or  replacement. 

To  insure  that  system  accuracy  will  stay  within  the  requirements6,  the  updates  of  the  navigation  parameters 
(orbit  elements  and  clocks)  must  occur  within  the  allowed  limits  of  the  degradation  of  the  parameters 
themselves.  These,  in  turn,  are  related  to  the  intrinsic  orbit  stability  (this  takes  into  account  gravitational 
field  asymmetries,  luni-solar  perturbations,  planetary  perturbations,  solar  pressure  effects,  etc.)  and  to  the 
clock  stability.  An  example  of  orbit  and  clock  degradation  as  they  appear  to  the  user  is  shown  in  Table  1 
and  in  Figure  6  (from  [4]).  Since  the  overall  accuracy  is  the  sum  of  these  two  causes  of  degradation  (orbit 
perturbations  plus  clock  stability),  a  stable  clock  will  help  to  maximize  the  update  period  for  a  given 
requirement.  During  each  update  interval,  measurements  are  collected  to  generate  a  new  set  of  parameters 
to  be  uploaded  at  the  end  of  the  data  collection  and  processing  sequence.  The  parameters  degradation  rate 
allows  the  setting  of  a  maximum  allowable  update  time  interval  between  updates.  Two  approaches  can  be 
selected. 

To  maintain  the  navigation  accuracy  of  the  system,  orbit  and  time  parameters  need  to  be  periodically 
updated,  since  both  tend  to  drift  with  time  and  to  increase  the  related  errors.  The  periodical  maintenance  is 
based  on  frequent  updates  of  the  navigation  parameters  and  less  frequent  tuning  of  system  parameters  or 
hardware  calibration  or  replacement. 

In  principle,  as  far  as  the  clocks  are  concerned,  the  basic  and  simpler  approach  will  be  to  select  an  update 
time  interval  such  that  the  clock  contribution  to  the  total  system  error  will  be  within  specified  limits  and, 
therefore,  negligible.  This  approach  has  the  advantage  that  modeling  of  the  clock  is  not  necessary. 
Enough  measurements  are  taken  during  each  measurement  interval  (less  or  equal  to  the  interval  between 
updates)  to  solve  completely  for  the  orbit  parameters  and  clock  offsets.  The  latter  are  considered  constant 
over  the  measurement  and  prediction  intervals,  since  the  clock  stability  guarantees  that  the  contribution  of 
the  clock  drift  and  instability  will  be  negligible  (only  the  frequency  offset  needs  to  be  estimated). 


6  Requirements  here  and  in  the  following  are  intended  for  the  “normal”  operation  of  the  system;  survivability 
considerations  when  normal  operations  are  disrupted  will  pose  additional  burdens  on  clocks  long-term  stability,  but 
the  relevant  considerations  are  limited  to  military  systems  in  particular  phases  of  their  operational  life. 
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Table  1  -  Projected  Block  IIF  Performance  in 
Crosslink  Navigation  Update  Mode  (Age  of  Data  of 
3  hours);  from  [4], 
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Figure  4  -  Predicted  Block  IIF  PPS  URL  vs.  Age  of 
Data  in  Crosslink  Navigation  Update  Mode  [4j 


The  possible  drawback  of  this  method  follows  from  the  consideration  that  measurements  taken  in 
successive  update  intervals  are  independent  except  lor  the  orbit,  which  is  propagated.  Therefore,  residual 
errors  may  accumulate  in  the  clocks  terms  without  bounds.  This  may  result  in  a  long-term  orbit  prediction 
which  is  very  uniform  but  not  very  accurate,  since  any  drill  in  the  orbital  parameters  may  simply 
accumulate  in  the  clock  terms,  which  are  considered  independent  from  one  measurement  period  to  the 
successive  one. 

Given  the  fact  that  atomic  clocks  are  quite  stable  [5j  and  their  behavior  can  be  predicted  over  long  time 
spans,  the  observations  during  one  measurement  period  may  be  complemented  by  the  prediction  of  the 
clock  behavior  over  the  past.  The  result  is  that  the  overall  system  errors  are  now  bound  by  the  clock 
stability  over  long  time  intervals,  observations  are  supplemented  by  predictions  which  add  additional 
information  and  help  to  detect  clock  malfunctions  in  the  system.  The  maintenance  operation  now  reverts 
to  a  classic  predictor-corrector  process,  where  the  prediction  is  provided  to  the  users  tor  navigation 
purposes,  but  is  also  used  by  the  system  for  the  correction  of  the  orbit  and  clock  parameters.  In  the  new 
situation,  the  clock  estimates  must  be  extended  beyond  Ehe  update  interval,  where  the  clock  contribution  to 
the  total  error  was  considered  negligible.  Then,  a  careful  modeling  of  the  clock(s)  behavior  is  mandatory, 
to  bound  the  errors  to  the  measurements  uncertainties  alone  and  not  to  introduce  additional  errors  due  to 
the  modeling  and  estimation  process.  Luckily,  the  systematic  and  stochastic  behavior  of  stable  oscillators 
has  been  intensively  studied  in  the  past  40  years,  so  current  models  are  at  least  as  good  as  the  clocks 
themselves. 


CLOCK  ERROR  MODEL 


A  good  clock  error  model  is  of  paramount  importance  in  the  design  and  operational  maintenance  phases  of 
a  satellite-based  navigation  system  For  our  purposes,  a  clock  is  a  numerical  integrator  of  the  periodic 
signal  produced  by  a  stable  oscillator;  there  lore,  the  timekeeping  function  is  implemented  by  counting 
cycles,  i.e.,  by  numerically  integrating  the  output  frequency  of  the  oscillator. 
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Frequency  accuracy  is  of  paramount  importance  for  an  accurate  clock,  while  stability  insures  that  accuracy 
is  maintained  over  time,  thereby  insuring  the  uniformity  and  reproducibility  of  the  time  scale. 

Stable  atomic  clocks  are  affected  by  systematic  and  stochastic  errors  that  must  be  accurately  modeled  and 
characterized  to  insure  optimum  performance  in  the  use  of  the  device  [I].  Generally,  the  basic  model  for  a 
stable  clock  includes  three  parameters: 

1)  an  initial  time  offset,  art=At(t0),  which  represents  the  inaccuracy  in  the  setting  of  the  initial  time  at 

2)  an  initial  frequency  offset,  a*,  which  represents  the  inaccuracy  in  selling  the  frequency  of  the  clock  at 
the  time  this  parameter  is  customarily  represented  as  normalized  with  respect  to  the  nominal 
frequency  of  the  oscillator  f0,  and  in  this  case  it  is  referred  to  as  the  initial  fractional  frequency  offset  of 
the  oscillator,  A/(/0)//0  ; 

3)  a  frequency  drift  a2  (a  parameter  referred  to  as  ageing  in  crystal  oscillators),  which  represents  the 
change  in  time  of  the  frequency  offset.  This  parameter  is  generally  assumed  constant  (or  slowly 
varying  with  time)  for  all  practical  purposes. 

Considering  the  clock  as  a  numerical  integrator,  the  above  three  terms  can  be  regarded  as  the  integration 
constants  for  the  cascaded  integrations  that  made  up  the  oscillators  and  the  clock:  these  are  schematically 
shown  in  Figure  5, 


Here  we  have  separated  the 
oscillator  model  and  the 
clock  model,  each  with  its 
own  integrators.  The  accum- 
ulated  time  error  Ai  (t)  al  the 
time  t  is  a  function  of  the 
above  parameters  and  the 
integration  functions,  and  can 
be  expressed  by  the  well- 
known  formula: 


Figure  5  -  Clock  systematic s  model,  including  environmental  and  random 
effects 


r  f  t* 


(3) 


Vo 


At(l)  =  au  +  Jo,  +  J  ja2  ■  dr'  ■  dt"  =  a0  +  at  ■  (t  ~t0 )  +  —  •  a2  ■  {t-tQ  )2 


This  model,  so  far.  includes  only  the  systematic  terms  intrinsic  to  the  oscillator  and  clock,  but  not  external 
systematic  effects  that  affect  any  additional  departure  from  the  nominal  frequency  and  phase  settings. 
These  environmental  effects  affect  the  frequency  of  the  oscillator,  through  temperature  dependence. 
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acceleration,  and  magnetic  field  sensitivity,  or  voltage  and  load  frequency  pulling  effects.  Moreover,  the 
clock  integration  may  be  affected  too,  for  instance  by  delay  changes  in  the  electronics  induced  by 
temperature  effects.  Adding  these  additional  environmental  effects,  the  resulting  model7  is  graphically 
shown  in  the  block  diagram  (Figure  5). 

The  noise  processes  affecting  the  output  frequency  of  the  oscillator  dominate  the  stochastic  behavior  of  the 
clock.  These  processes  have  been  intensively  studied  during  the  past  40  years,  and  are  characterized  by  a 
measurement  of  frequency  instability  expressed  as  spectral  density  of  the  phase  noise  (in  the  frequency 
domain)  and  Allan  variance  (in  the  time  domain)[6].  Since  for  a  clock  we  are  usually  interested  in  the  time 
domain  aspects  of  frequency  instabilities,  equation  (3)  is  completed  to  account  for  stochastic  effects  as 
follows  (including  environmental  effects): 

A t(t)  =  a0  +al  '(t-t0)+^a2\t-t0)2  + 

<0  k  J 

where  the  upper  line  in  (4)  accounts  for  the  systematics,  the  first  term  of  the  second  line  accounts8  for  the 
nonstationary  noise  processes  of  the  oscillator9,  the  second  term  for  the  environmental  effects  affecting  the 
frequency  (the  k  parameters  include  temperature  sensitivity,  magnetic  sensitivity,  load-pulling  effects, 
relativity,  etc.),  and  the  last  term  the  delay  (phase)  fluctuations  induced  by  temperature  in  the  clock 
(counter’s)  electronics.  More  details  on  the  clock  model  can  be  found  in  [3,6-10]. 

In  practice,  two  steps  are  required  to  maintain  a  timing  system  when  the  data  are  acquired  in  batches  and 
the  clock  state  vector  is  estimated  at  regular  intervals  (as  in  the  case  of  navigation  systems): 

1.  characterization  of  the  clocks  behavior,  i.e.,  estimation  via  measurements  of  the  relevant  parameters 
entering  Eq.  (4),  namely: 

time  offset  a<>  at  some  (initial)  time  to  -  for  all  purposes;  to  can  be  selected  as  completely  arbitrary; 
fractional  frequency  offset  ai  at  to, 

frequency  drift  a2,  assuming  that  this  term  is  constant  to  avoid  a  further  integration  in  the  model10; 

frequency  stability,  expressed  as  the  Allan  variance  <72  (r)over  the  time  interval  of  interest  (in  this 
case  the  time  required  to  carry  on  the  measurements  used  to  characterize  the  device) 
sensitivity  parameters  pk  of  the  clock  to  the  various  environmental  factors. 

2.  prediction  of  the  clock  behavior,  i.e.  forecasting  of  the  clock  error  A t{t)  in  time  over  some  time 


7  We  have  not  considered  random-walk  noise  processes,  since  they  generally  affect  the  oscillator  at  longer  integration 
times  than  those  considered  for  updating  in  the  normal  operations  of  a  navigation  system.  These  can  be  added  to 
the  model  if  required  for  simulation  purposes. 

8  This  is  not  a  very  precise  model  of  the  time  deviation  of  the  clock  due  to  the  frequency  stability  effects,  and  better 
predictors  can  be  used  when  the  various  noise  processes  contributions  are  taken  into  account  separately. 

9  V(x)  and  ti(t)  are  the  contributions  due  to  the  white  frequency  and  flicker  frequency  noise  processes  respectively. 

10  One-way  pseudorange  measurements  may  be  supplemented  by  two-way  data  for  the  estimation  of  ai  (and  a2  when 
needed)  for  the  MS  clocks  -  see  note  4  above. 
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interval  in  the  future,  based  again  on  the  clock  model  expressed  by  (4),  where  now  the  relevant 
parameters  are  based  on  the  estimate  { I )  above. 


By  considering  that  the  stochastic  components  are  given  by  gaussian  noise  with  zero  average,  the 
prediction  of  the  clock  error  A/(f)  comes  straightforward  from  Fq.  (4)  above.  The  uncertainty  of  the 
prediction  has  to  take  into  account  the  uncertainty  in  the  estimation  of  the  systematic  parameters,  as  well  as 
the  possible  covariances  between  estimates,  and  also  the  fact  that  the  estimation  of  the  stochastic 
components  is  influenced  and  influences  the  estimation  of  the  systematic  elements.  The  optimal  estimation 
of  these  components  is  at  the  core  of  the  current  debate  [9,10 1  Let's  concentrate  at  the  moment  only  on 
the  prediction  of  the  stochastic  components.  The  most  common  clock  noises  are  gaussian  and  the 

uncertainty  in  the  clock  error  A/(f)  prediction  may  be  related  to  the  Allan  deviation  by  appropriate 
relationships  [6,7,8],  In  case  of  white  and  flicker  EM,  for  example,  the  uncertainty  of  the  A/(/).  here 
named  by  assuming  may  be  written  as: 


(5) 


FFM  -  v 

=  o 


-/in  2 


where  a/'  w  (t)  stands  for  the  Allan  deviation  of  the  white  FM  at  time  /,  while  GyFFM  (t)  stands  for  the 

flicker  FM,  The  use  of  these  relationships  will  help  in  the  following  section,  where  the  prediction 
uncertainty  of  different  clocks  and  the  synchronization  update  rate  will  be  discussed. 


The  estimation  and  prediction  phases  will  generally  overlap,  since  the  clock  error  is  propagated  during  the 
acquisition  of  the  batch  of  measurements  required  to  characterize  the  clock  and  enable  the  time  error 
estimation  over  the  next  prediction  interval.  At  the  epoch  in  which  a  new  estimate  is  generated,  the  clock 
error  al  the  same  epoch  can  be  compared  with  the  clock  error  previously  estimated  and  the  inaccuracy  of 
the  characterizat  ion/predict  ion  process  can  be  assessed,  to  validate  both  the  measurement  process  and  the 
stale  model. 


Estimation/prediction  methods  generally  fall  in  three  main  classes: 

1.  ARMA  (Auto-Regressive  Moving  Average)  /ARIMA  (Auto-Regressive  Integrated  Moving 
Average)  methods  [2] 

2.  (Fixed)  filters  methods,  based  on  some  kind  of  functional  fitting,  generally  polynomial  [2,12  j 

3.  (Adaptive)  filters  methods,  based  on  Kalman  filters. 
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Irrespective  of  the  method,  the 
process  requires  a  data 
collection  interval*  wrhere 
measurements  are  collected  as 
a  time-ordered  series,  to 
generate  a  set  of  prediction 
parameters  that  constitute  l he 
system  state  description  over 
the  next  prediction  interval. 

The  frequency  stability  of  the  oscillator  plays  a  key  role  not  only  in  the  prediction  of  the  clock  behavior, 
but  also  in  the  estimation  of  the  systematic  parameters  at  the  end  of  each  measurement  interval.  As  a 
matter  of  fact,  the  frequency  stability  limits  the  estimation  of  the  systematic  parameters  of  the  clock,  and 
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the  result  is  that  the  predicted  A/  #  0  at  the  beginning  of  the  prediction  interval.  Since  the  aim  is  to 
minimize  A/  at  the  end  of  the  prediction  interval  when  a  new  batch  of  data  has  been  acquired  and  a  new 
dock  state  estimates  is  available,  two  constraints  should  be  met  to  minimize  A t : 

•  the  length  of  the  measurements/prediction  periods  should  be  made  as  short  as  possible,  taking  into 
account  all  operational  constraints,  since  during  the  prediction  the  error  increases  with  time 

*  the  length  of  the  measurements/prediction  periods  should  be  made  long  enough  to  allow  the  estimation 
of  the  clock  parameters  with  sufficient  accuracy. 


NEW  OSCILLATORS 

Clearly,  a  tradeoff  of  the  above  conditions  will 
highlight  the  necessity  to  reach  the  flicker 
frequency  floor  of  the  oscillator  as  soon  as 
possible,  to  allow  a  precise  characterization  of 
the  clock  systematic  parameters  in  the  shortest 
possible  time,  again  assuming  that  operational 
constraints  will  not  pose  further  limitations. 

Further,  for  shorter  measurement/prediclion 
intervals,  the  clock  model  can  be  simplified  by 
dropping  terms  that  have  negligible  effects  over 
the  measurement/prediction  interval. 


Figure  6  -  Frequency  instability  for  various  atomic 
oscillators,  including  the  proposed  CPT  Rb  maser 


However,  the  possibility  that,  for 
whatever  operational  reason,  the 
prediction  interval  may  be  extended 
during  contingency  situations,  may 
require  careful  considerations  in  the 
definition  of  a  simplified  model  (and 
also  in  the  choice  of  oscillators  with 
no  performance  margin). 

H-masers,  both  of  the  active  and 
passive  types,  are  characterized  by  an 
exceptional  short-  to  medium-term 
frequency  stability,  and  are  considered 
in  the  baseline  for  spaeebome 
oscillators  in  newr  systems,  like 
Galileo  (Figure  6), 
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However,  ihey  are  complex  devices  and  require  vacuum  pumps  and  hydrogen  reservoirs  LhaL  may  limit 
long-term  reliability,  even  if  outstanding  results  have  been  presented  in  recent  times.  An  interesting  device 
has  been  recently  proposed:  the  coherent  population  trapping  (CPT)  maser  [11].  To  maintain  the 
navigation  accuracy  of  the  system,  orbit  and  lime  parameters  need  to  be  periodically  updated,  since  both 
tend  to  drift  with  lime  and  to  increase  the  related  errors.  The  periodical  maintenance  is  based  on  frequent 
updates  of  the  navigation  parameters  and  less  frequent  tuning  of  system  parameters  or  hardware  calibration 
or  replacement.  This  oscillator  is  based  on  a  low  density  vapor  of  Cs  or  Rb  contained  in  a  cell  with  buffer 
gas  and  coupled  to  a  microwave  cavity  tuned  to  the  atomic  ground  state  hyper  fine  frequency. 

The  atomic  sample  is  excited  by  two  laser  fields  that  couple  the  two  ground  state  hyper  fine  levels  to  the 
excited  PW2  or  P m  states. 


This  A  scheme  excitation  process 
leads  to  the  presence  of  a  dark  line  in 
the  fluorescence  spectrum  and  to  an 
oscillating  magnetization  responsible 
for  a  coherent  microwave  emission  at 
the  hyperfme  frequency.  Both  signals 
may  be  used  as  atomic  reference  as 
shown  in  the  experimental  setup 
shown  in  Figure  7.  The  laser  carrier  is 
phase-modulated  by  an  electro-optic 
modulator  and  the  two  first  sidebands 
provide  the  radiation  fields  required 
for  the  excitation  of  the  atomic 
ensemble.  Figure  8  -  Time  prediction  uncertainty  for  various  atomic 

oscillators,  including  the  proposed  CPT  Rb  maser 


Time  prediction  uncertainty  (1  sigma) 


The  high  degree  of  symmetry  of  this  excitation  scheme  strongly  reduces  the  light-shift  effect  that  is 
considered  the  main  source  of  frequency  instability  for  the  optically  pumped  passive  standards  in  the  short- 
lo-medium  term. 

Moreover,  the  microwave  signal  is  observed  without  any  background  and  allows  a  strong  reduction  of  the 
laser  amplitude  fluctuations  transferred  on  the  output  signal.  The  timing  performance  of  this  device, 
resulting  from  the  projected  frequency  stability  (Figure  6),  stand  midway  between  the  classic  optically 
pumped  Rb  passive  cell  and  the  passive  H -maser,  providing  a  reasonable  option  for  future  spaceborne 
clocks. 
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